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We discuss properties of heavy quarkonium states at high temperatures based on lattice 
QCD and potential models. We review recent progress made in lattice calculations of 
spatial static quark anti-quark correlators as well as quarkonium correlators in Euclidean 
time. Recent developments in effective field theory approach and potential models are 
also discussed. 



1. Introduction 

There was considerable interest in the properties and the fate of heavy quarkonium 
states at finite temperature since the famous conjecture by Matsui and Satz [1]. It 
has been argued that color screening in medium will lead to quarkonium dissocia- 
tion above deconfinement, which in turn can signal quark gluon plasma formation 
in heavy ion collisions. The basic assumption behind the conjecture by Matsui and 
Satz was the fact that medium effects can be understood in terms of a temperature 
dependent heavy quark potential. Color screening makes the potential exponentially 
suppressed at distances larger than the Debye radius and it therefore cannot bind 
the heavy quark and anti-quark once the temperature is sufficiently high. Based 
on this idea potential models at finite temperature with different temperature de- 
pendent potentials have been used over the last two decades to study quarkonium 
properties at finite temperature (see Ref. [2] for a recent review). It was not until 
recently that effective field theory approach, the so-called thermal pNRQCD, has 
been developed to justify the use of potential models at finite temperature [3]. This 
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approach, however, is based on the weak couphng techniques and wiU be discussed 
in the next section. To understand the non-perturbative aspects of color screening 
lattice calculations of the spatial correlation functions of static quarks are needed. 
Recently a lot of progress has been made in this direction which will be the topic of 
section |4l To prepare the reader for this in section [3] we review the basics of lattice 
gauge theory. 

In principle it is possible to study the problem of quarkonium dissolution without 
any use of potential models. In-medium properties of different quarkonium states 
and/or their dissolution are encoded in spectral functions. Spectral functions are 
related to Euclidean meson correlation functions which can be calculated on the 
lattice. Reconstruction of the spectral functions from the lattice meson correlators 
turns out to be very difficult, and despite several attempts its outcome still remains 
inconclusive. One remarkable feature of the studies of the lattice meson correlators 
is their small temperature dependence despite the expected color screening. This 
seems to be puzzling. We will discuss the possible resolution of this puzzle in section 
m while the current status of the lattice calculations of the Euclidean correlators 
and the corresponding meson spectral functions will be presented in section [5l The 
summary and outlook will be given in section [7l 

2. pNRQCD at finite temperature 

There are different scales in the heavy quark bound state problem related to the 
heavy quark mass m, the inverse size ~ mv and the binding energy mv'^. Here v is 
the typical heavy quark velocity in the bound state and is considered to be a small 
parameter. Therefore it is possible to derive a sequence of effective field theories 
using this separation of scales (see Refs. [4, 5] for recent reviews). Integrating out 
modes at the highest energy scale ~ m leads to an effective field theory called non- 
relativistic QCD or NRQCD, where the pair creation of heavy quarks is suppressed 
by powers of the inverse mass and the heavy quarks are described by non-relativistic 
Pauli spinors [6]. At the next step, when the large scale related to the inverse size is 
integrated out, the potential NRQCD or pNRQCD appears. In this effective theory 
the dynamical fields include the singlet S(r, R) and octet 0(r, R) fields correspond- 
ing to the heavy quark anti-quark pair in singlet and octet states respectively, as 
well as light quarks and gluon fields at the lowest scale ~ mv^. The Lagrangian of 
this effective field theory has the form 




(1) 
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Here the dots correspond to terms which are higher order in the muhipole expansion 
[5] . The relative distance r between the heavy quark and anti-quark plays a role of a 
label, the light quark and gluon fields depend only on the center-of-mass coordinate 
R. The singlet Vs{r) and octet Vo{r) heavy quark potentials appear as matching 
coefficients in the Lagrangian of the effective field theory and therefore can be 
rigorously defined in QCD at any order of the perturbative expansion. At leading 
order 

and Va = Vb = I- 

The free field equation for the singlet field is 



2 



ido + ^-Vs{r) 
m 



S(r,R)=0, (3) 



i.e. has the form of a Schrodinger equation with the potential Vs{r). In this sense, 
potential models emerge from the pNRQCD. Note, however, that pNRQCD also 
accounts for interaction of the soft gluons which cannot be included in potential 
models, i.e. it can describe retardation effects. 

One can generalize this approach to finite temperature. However, the presence of 
additional scales makes the analysis more complicated [3] . The effective Lagrangian 
will have the same form as above, but the matching coefficients may be temperature 
dependent. In the weak coupling regime there are three different thermal scales : T, 
gT and g^T. The calculations of the matching coefficients depend on the relation 
of these thermal scales to the heavy quark bound state scales [3]. To simplify the 
analysis the static approximation has been used, in which case the scale mv is 
replaced by the inverse distance 1/r between the static quark and anti-quark. The 
binding energy in the static limit becomes Vo — Vs — Nag/ (2r). When the binding 
energy is larger than the temperature the derivation of pNRQCD proceeds in the 
same way as at zero temperature and there is no medium modifications of the heavy 
quark potential [3]. But bound state properties will be affected by the medium 
through interactions with ultra-soft gluons, in particular, the binding energy will 
be reduced and a finite thermal width will appear due to medium induced singlet- 
octet transitions arising from the dipole interactions in the pNRQCD Lagrangian [3] 
(c.f. Eq. ([1])). When the binding energy is smaller than one of the thermal scales 
the singlet and octet potential will be temperature dependent and will acquire 
an imaginary part [3]. The imaginary part of the potential arises because of the 
singlet-octet transitions induced by the dipole vertex as well as due to the Landau 
damping in the plasma, i.e. scattering of the gluons with space-like momentum off 
the thermal excitations in the plasma. In general, the thermal corrections to the 
potential go like (rT)" and (mur)" [3], where mo denotes the Debye mass. Only 
for distances r > 1 /m^i there is an exponential screening. In this region the singlet 
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potential has a simple form 

Vs{r) = -Cf ^ e-"^-^ + zCf a.T ^ T dx ^^"("^^y^) _ Cp a, {mo + iT) , 
r rm£i Jq (x^ + 1)"^ 

Cf = (N^ - 1)/(27V) (4) 

The real part of the singlet potential coincides with the leading order result of 
the so-called singlet free energy [7] . The imaginary part of the singlet potential in 
this limit has been first calculated in [8]. For small distances the imaginary part 
vanishes, while at large distances it is twice the damping rate of a heavy quark [9]. 
This fact was first noted in Ref. [10] for thermal QED. 

The effective field theory at finite temperature has been derived in the weak cou- 
pling regime assuming the separation of different thermal scales as well as Aqcd- 
In practice the separation of these scales is not evident and one needs lattice tech- 
niques to test the approach. Therefore section [3] will be dedicated to the study of 
static quarks at finite temperature on the lattice. To prepare the reader for this 
discussion some basics of the lattice gauge theory will be given in the next section. 



3. Basics of lattice gauge theory 

To study non-perturbative aspects of QCD we use lattice gauge theory [11]. In this 
formalism a field theory is defined in a gauge- invariant way on a discrete space-time 
domain. This serves at least two purposes: a) to provide an ultra-violet cut-off for 
the theory, restricting highest momentum to n/a (a being the lattice spacing), and 
b) to evaluate the path integrals in the Euclidean formulation stochastically using 
importance sampling. 

On the lattice the fundamental degrees of freedom of a theory with local SU (TV) 
gauge symmetry are fermion fields ipx that reside on the sites of the lattice and carry 
flavor, color and Dirac indeces, which we suppress through the most of this paper, 
and gauge, bosonic degrees of freedom that in the form of SU{N) matrices Ux,f_i 
reside on links. Sites on a four-dimcnsional lattice are labeled with x = {x,t). 

The theory is defined by the partition function 

Z = J DUDtpD^expi-S) (5) 

where the action 

S=Sg + Sf (6) 

contains gauge, Sg and fermionic, Sf parts. The latter part is bi-linear in fields and 
has the form 

Sf = {pM^ (7) 

where M is the fermion matrix. In the simplest formulation the lattice gauge action 
can be written as 



^.-/?E(l-^TrC/p), 



(8) 
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where Up is the so-caUed plaquette, a product of hnk variables along the elementary 
square and (3 = 'iN/g'^ with being the bare gauge coupling. This is the Wilson 
gauge action [11]. The explicit form of the fermion action that is often used in 
lattice QCD calculations will be discussed in section [5?2l 
The expectation value of an operator O is given then by 

{O}^ ^ J DUD'4)Di:dexp{~S). (9) 

Integration over the fermion fields (which are Grasmann variables) can be carried 
out explicitly: 

Z = J DUdetM[U]exp{-Sg) = J DU cxp{-Seff ), (10) 

where Seff — Sg — \ndet M[U] is the effective action. The fermion determinant 
det M[U] describes the vacuum polarization effects due to the dynamical quarks 
and makes the effective action non-local in gauge variables. For this reason sim- 
ulations with dynamical quarks are very resource demanding and the quenched 
approximation is often employed, where detM[J7] is set to 1. 

To evaluate the path integral ([9]) stochastically, an ensemble of Nu gauge config- 
urations, weighted with cxp{—Seff), is generated using Monte Carlo or Molecular 
Dynamics techniques. The expectation value of the operator is then approximated 
by the ensemble average: 

, Nu 

(0)^J^T.O^iU), (11) 
i—1 

where Oi{U) is the value of the operator O, calculated on i-th configuration. (When 
the operator O depends explicitly on the quark fields extra factors of appear 
as shown for a meson correlator below.) 

Consider a meson (quark anti-quark pair) operator of a general form 

J{x,y;t) = i^{x,t)m{x,y;t)ij{y,t), (12) 

where V determines the spin structure and hi is a, gauge connection that corresponds 
to the excitations of the gluonic field. Dirac and color indcccs in are suppressed. 

The propagation of such meson from time t = to t is described by the corre- 
lation function 

(J(a-i,yi;0)J(x2,y2;i)> = ^ j DUD^pD^expi-S) 

X i^(xi,0)rZ^(xi,yi;0)^(yi, 0)^(^2, i)rtZ^^(x2,y2;0V'(^2,i). (13) 
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Again, integration over the quark fields can be carried out resulting in 

{J{xi,yi;0)J{x2,y2;t)) = 

(Tr [M-\x2,t; xi,0)TU{xi,yi;0)Ar\yi,0; ys, t)r^U^X2,y2; t)] ) 
- (Tr [M-\y,,0; fi, 0)m(xi, yi; 0)]) 

X {TT[M-\x2,t;y2,t)T^UHx2,y2]t)]). (14) 

The inverse of the fermion matrix, M^^ has meaning of a fermion propagator. 

Gauge transporters U{x,y;t) can be taken as weighted sums of diflFerent paths 
connecting points x and y. By choosing paths of certain shape or combinations of 
different paths it is possible to achieve a better overlap of the meson operator with 
a given state. One of the possibilities is to construct the gauge transporters by using 
APE smearing [12] on spatial links: a link variable Ux.fj. is replaced by a weighted 
average of itself and a sum of the 3-link paths connecting the same sites as Ux.^' 

U.,^ ^ K..^ = (1 - 6c)C/.,^ + c J2 Ux,,Ux+,.,,Ul^.^,^. (15) 

This procedure can be applied iteratively. Then a gauge transporter U{x,y;t), taken 
as a product of smeared links, is equivalent to a weighted sum of differently shaped 
paths connecting the sites x and y. 

In the following we will consider meson correlators at finite temperature. The 
finite temperature is introduced by compactifying the the Eucliden time direction, 
i.e. T — l/{Nra) with Nr being the number of temporal time slices. Gauge fields and 
fermion fields obey periodic and anti-periodic boundary conditions in the temporal 
direction. In the next section we consider a spinless static quark anti-quark pair, 
r = I that can propagate only in time. In this case second term in ([H]) vanishes. 
In Sec. [5] we consider local meson operators with Xi = yi, i ~ 1,2 which means 
U{x,y\t) = I. 

4. Correlation functions of static quarks in lattice gauge theory 
4.1. Static meson correlators 

Consider static (infinitely heavy) quarks. The position of heavy quark anti-quark 
pair is fixed in space and propagation happens only along the time direction. In this 
limit the second term on the right hand side of (fT4|) vanishes. We are interested in a 
spinless state and set F = I in this section. With respect to the color the meson can 
be in a singlet or adjoint state. These states are described by the following gauge 
connections 

U{x,y;t)^U{x,y;t), (16) 
W{x,y;t) = C/(x,xo;t)T"[/(xo,y;t), (17) 

where U{x, y) is a spatial gauge transporter, - the product of the gauge variables 
along the path connecting x and y, xq is the coordinate of the center of mass of the 
meson and T" are the SU (TV) group generators. 
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The meson operators are given then by Eq. with xi = X2, yi — 2/2 for static 
quarks 

J{x, y: t) = i,{x, t)U{x, y; t)i^{y, t), (18) 
r{x,y;t) = 4,{x,t)U{x,xo;t)T''U{xo,y;t)ib{y,t). (19) 

Substituting expressions (|18p into eq. p4p and noting that for a static quark 
the propagator M~^{x,0; x,t) ~ L{x), where the temporal Wilson line L{x) ~ 
Ylt^Q ^ t^(x,t),o with C/(x,t),o being the temporal links, we get for the meson correla- 
tors at i = 1/T: 

Gi(r,T) = l(J(f,z7;0)J(f,y;l/T)) 

= l(Tr [LHx)U{x,y;0)L{y)U\x,y,l/T)]), (20) 

Ga{r,T) ^ {J''{^:y;0)r{S;,y;l/T)) 

a=l 

= ^^2^(Tr£t(x)TrL(y)) 

- ^(^2 _ 1) [£^(x)^(x, y; 0)L{y)uHx,y, 1/T)]), (21) 

The correlators depend on the choice of the spatial transporters U (x, y; t). Typ- 
ically, a straight line connecting points x and y is used as a path in the gauge trans- 
porters, i.e. one deals with time-like rectangular cyclic Wilson loops. This object 
has been calculated at finite temperature in hard thermal loop (HTL) petrurbation 
theory in context of perturbative calculations of the singlet potential introduced in 
the previous section and quarkonium spectral functions [8, 13, 14]. In the special 
gauge, where U{x, y;t) = 1 the above correlators give the standard definition of the 
singlet and adjoint free energies of a static QQ pair 

exp(-^^i(r,r)/T) - l(Tr[Lt(a:)L(y)]), (22) 

exp(-F,(r,T)/T) = -^{TTL\x)TTL{y)) 

1 



-^-^^iTr[LHx)Liy)]). (23) 



The singlet and adjoint free energies can be calculated at high temperature in 
leading order HTL approximation [7] resulting in 

Fi{r,T) = exp(-mz,r) — , (24) 

Pair^T) . ^^^exp(-™..) - (25) 
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with niD = gT^y{N/3 + Nf/6) being the leading order Debye mass and Nf is the 
number of quark flavors. At this order Fi and Fa are gauge independent or, in other 
words, do not depend on the choice of the parallel transporters U{x, y; t). Note that 
at small distances {rmo ^1) the singlet free energy 

^^i(r,T)^-^^^ (26) 

is temperature independent and coincides with the zero temperature potential, 
while the adjoint free energy 

Fa{r,T) ^ - ^asmn (27) 

depends on the temperature. 

The physical free energy of a static quark anti-quark pair, i.e. the one related 
to the work that has to be done to separate the static charges by certain distance 
is given by the thermal average of the singlet and adjoint free energies [15] 

1 N'^ -1 

exp(-F(r,T)/T) = _ exp(-i^i(r, T)/T) + exp(-i^,(r, r)/r) 

= ^(Tr[L(x)TrL(y)]) ^ ^G{r,T). (28) 

This quantity is explicitly gauge independent. In leading order HTL approximation 
the free energy is 

(N^ - 1) 

nr.T) ^ -i^^-^exp(-2™^r). (29) 

The 1/r^ behavior is due to partial cancellation between the singlet and adjoint con- 
tribution [15, 16] and has been confirmed by lattice calculations in the intermediate 
distance regime above the deconfinemcnt transition [17, 18]. 

Using the transfer matrix one can show that in the confined phase 



Gi(r,T) = ^c„(r)e-^"(-^)/^ (30) 
G(r,T) = f]e-^"('-'^)/^, (31) 

n=l 

where En are the energy levels of static quark and anti-quark pair [19]. The coef- 
ficients c„(r) depend on the choice of U{x, y; t) entering the static meson operator 
in Eqs. (|18m9p . Since the color averaged correlator G{r,T) corresponds to a gauge 
invariant measurable quantity it does not contain c„. The lowest energy level is the 
usual static quark anti-quark potential, while the higher energy levels correspond 
to hybrid potentials [20-23]. Using multi-pole expansion in pNRQCD one can show 
that at short distances the hybrid potential corresponds to the adjoint potential up 
to non-perturbative constants [24] . Indeed, lattice calculations of the hybrid poten- 
tials indicate a repulsive short distance part [20-23]. Furthermore, the gap between 
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the static potential and the first hybrid potential can be estimated fairly well at 
short distances in perturbation theory [25]. 

If ci = 1 the dominant contribution to Ga would be the first excited state E2, 
i.e. the lowest hybrid potential which at short distances is related to the adjoint 
potential. In this sense Ga is related to static mesons with quark anti-quark in 
adjoint state. Numerical calculations show, however, that ci is r-dependent and in 
general ci(r) 7^ 1. Thus Ga also receives contribution from Ei [19]. The lattice data 
suggests that ci approaches unity at short distances [19] in accord with expectations 
based on perturbation theory, where ci = 1 up to 0{a^) corrections [24]. Therefore 
at short distances, r ^ 1/T the color singlet and color averaged free energy are 
related F(r,T) = Fi{r,T) + Tlii{N^ - 1). 

In the following we consider SU (2) and SU (3) gauge theories and refer to the 
adjoint state as triplet and octet, correspondingly. 

4.2. Lattice results on static meson correlators 

Correlation function of static quarks have been extensively studied on the lattice 
since the pioneering work by Mclerran and Svetistky [15]. Most of these studies, 
however, considered only the color averaged correlator, i.e. the correlation function 
of two Polyakov loops (for the most complete analysis see Ref. [26] and references 
therein). Since both color singlet and color octet degrees of freedom contribute to 
the Polyakov loop correlator it has large temperature dependence even at short 
distances and it is not a very useful quantity if we want to learn something about 
quarkonium properties at high temperatures. As this has been pointed out in Refs. 
[27, 28] we need to know the quark anti-quark interactions in the singlet channel in 
order to learn about in- medium quarkonium properties. Therefore in recent years 
the singlet correlator has been computed on the lattice in SU{2) and SU{3) gauge 
theories [18, 29-32] as well as in full QCD with 3 and 2 flavors of dynamical quarks 
[33, 34]. Preliminary results also exist for 2-1-1 flavor QCD with physical value of 
the strange quark mass and light m, d-quark masses corresponding to pion mass of 
about 220MeV [35, 36]. All these calculations use Coulomb gauge definition of the 
singlet correlator, i.e. the definition ((22|) with Coulomb gauge fixing. 

The numerical results for SU(3) gauge theory are presented in Figure [1] and 
compared with the zero temperature potential. Here we used the string Ansatz for 
the zero temperature potential 



since this form gives a very good description of the lattice data in SU{3) gauge 
theory [37] . To convert lattice units to physical units the value y/a = 420MeV has 
been used for the string tension. To remove the additive renormalization in the 
singlet free energy the lattice data have been normalized to V{r) given by Eq. ([32]) 
at the shortest distance available. As one can see from Figure [1] the singlet free 
energy is temperature independent at short distances and coincides with the zero 



V{r) 



TT 



(32) 
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0.5 1 1.5 2 2.5 3 

Fig. 1. The color singlet free energy singlet free energy in quenched QCD {SU(3) gauge theory) 
calculated in Coulomb gauge [29, 31, 32]. The solid line shows the zero temperature potential. 

temperature potential. Below the deconfinement transition, T < Tc, it rises linearly, 
indicating confinement. The string tension at finite temperature is smaller than the 
zero temperature string tension a. At large distances the singlet free energy is the 
same as the free energy determined from the Polyakov loop correlators and therefore 
the finite temperature string tension agrees with findings of Ref. [26]. Medium 
effects set in at distance rmed — 0.4 im/{T/Tc) and we see exponential screening 
at distance r > 1/T [32]. In the intermediate distance regime 0.5 fm < r < 1.5 fm 
we see that the singlet free energy is enhanced relative to the zero temperature 
potential. This is not a real physical effect but an artifact of the calculations. Similar 
effect has also been seen in 5/7(2) gauge theory [18, 19, 30]. At these distances the 
color singlet correlator is sensitive to the value of the coefhcients c„(r) in Eq. ([30|) . 
Below we will show that the enhancement of the singlet free energy is due to the 
fact that ci(r) < 1 in this region. 

Numerical results for the singlet free energy in 2+1 flavor QCD are shown 
in Figure [2l In this case a different normalization procedure has been used. The 
additive renormalization has been determined at zero temperature for each value 
of the lattice spacing used in the finite temperature calculations by matching the 
zero temperature potential to the string Ansatz ([32|) at distance r = tq, with ro 
being the Sommer scale [38]. This additive renormalization then has been used 
for the singlet free energy. For detailed discussion of the calculation of the static 
potential and its renormalization see Ref. [39]. As in quenched QCD the singlet 
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0.5 1 1.5 2 2.5 



r[fm] 

Fig. 2. The color singlet free energy in 2+1 flavor QCD [35, 36] . The solid line is the parametrization 
of the lattice data on the zero temperature potential from Ref. [39]. 

free energy is temperature independent at short distances and coincides with the 
zero temperature potential. At distances larger than the inverse temperature it is 
exponentially screened. The novel feature of the singlet free energy in full QCD is 
the string breaking, i.e. the fact that it approaches a constant at large separation. 
This happens because if the energy in the string exceeds the binding energy of a 
heavy-light meson the static quark and anti-quark are screened due to pair creation 
from the vacuum. This happens at distances of about 0.8 fm according to the figure. 
As temperature increases the distance where the singlet free energy flattens out 
becomes smaller and for sufficiently high temperatures turns out to be inversely 
proportional to the temperature. Thus string breaking smoothly turns into color 
screening as temperature increases. 

4.3. Color singlet correlator in SU(2) gauge theory at low 
temperatures 

For a better understanding of the temperature dependence of the singlet correlator 
and its physical interpretation the values of the overlap factors c„(r) should be 
estimated. To extract the overlap factors we need to calculate the singlet correla- 
tors in the low temperature region, where only few energy levels contribute to the 
correlator. The calculations of the singlet correlator at low temperatures is difficult 
because of the rapidly deceasing signal to noise ratio. To overcome this difficulty it 
has been suggested to calculate Wilson loops with multilevel Luescher-Weisz algo- 
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1.5 
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Fig. 3. The color singlet free energy in SU{2) gauge theory below the deconfinement temperature 
at /3 = 2.5 calculated on 32^ X A^t lattices. Also shown is the T = potential. The inset shows 
the color singlet free energy from which the contribution from the matrix element T In ci has been 
subtracted. 



rithm [40] instead of the Coulomb gauge fixing [41]. In this case the color singlet 
correlator defined by Eq. (pO| is an expectation value of the gauge-invariant Wilson 
loop with gauge transporter U {x, y; t) being a product of the iteratively smeared 
spatial links. 

Numerical calculations have been performed in SU (2) gauge theory using stan- 
dard Wilson gauge action [41]. With the use of the multilevel algorithm it was 
possible to go down to temperatures as low as 0.32Tc not accessible in the previ- 
ous studies. It is well known that smearing increases the overlap with the ground 
state by removing the short distance fluctuations in the spatial links [42]. For this 
reason smearing also reduces the breaking of the rotational invariance to the level 
expected in the free theory. When no smearing is used the color singlet free energy, 
— T In Gi (r, T) shows a small but visible temperature dependence. The temperature 
dependence of the singlet free energy is significantly reduced when APE smearing 
is applied. The color singlet free energy for (3 = 2.5 and 10 APE smearings is shown 
in Fig. [3l As one can see from the figure the color singlet free energy shows very 
mild temperature dependence in the confined phase with noticeable temperature 
effects appearing only at T = O.QSTc- 

Consider the expansion (|30|) . The dominant contribution comes from the ground 
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state, so it is reasonable to fit the singlet correlator to the form 



Gi(r, r) = ci(r) exp(-£;i (r)/T). 



(33) 



This allows to extract the matrix element ci(r) using a simple exponential fit, which 
is shown in Fig. 21 When no APE smearing is used the value of Ci(r) strongly 
depends on the separation r. At small distances it shows a tendency of approaching 
unity as one would expect in perturbation theory. However, ci(r) decreases with 
increasing distance r. At large distance its value is around 0.3 — 0.5. Similar results 
for ci (r) have been obtained in the study of SU{2) gauge theory in 3 dimensions [19] . 
When APE smearing is applied the r-dependence of the amplitude ci (r) is largely 
reduced and its value is close to unity both for (3 — 2.5 and (3 = 2.7. For /? = 2.7 
we also see that increasing the number of smearing steps from 10 to 20 reduces the 
deviation of ci (r) from unity. 

As discussed in section 14.11 perturbation theory predicts that the deviation of 
Ci(r) from unity is of order a^. Therefore it can be made arbitrarily small by 
going to sufficiently small distances. It is known, however, that lattice perturbation 
theory converges very poorly. The main reason for this has been identified with 
the short distance fiuctuations of the link variables, which makes their mean value 
very different from unity [43]. Smearing removes these short distance fluctuations 
and this is the reason why ci (r) is much closer to unity when APE smearing is 
applied. Thus, almost the entire temperature dependence of the singlet free energy 
at distances 0.5 < ry/a < 2 is due to the deviation of ci from unity and can 
be largely reduced by applying APE smearing to the links in the spatial gauge 
connections. To further demonstrate this point in the inset of Figure [3] we show 
the results for Fi{r, T) +nnci(r). Clearly no temperature dependence can be seen 
in this quantity up to 0.95Tj,, where we see temperature dependence at distances 
r^fa > 1.5 corresponding to the expected drop of the effective string tension. 

4.4. Color singlet free energy in the deconfined phase 

The behavior of the color singlet free energy in the deconfined phase has been 
studied in Coulomb gauge [18, 29, 32-34] and from cyclic Wilson loops [41]. As 
discussed above at short distances it is temperature independent and coincides 
with the zero temperature potential. At large distances it approaches a constant 
Foo{T), which monotonically decreases with the temperature. The constant Foo{T) 
is the free energy of two isolated static quarks, or equivalently of a quark anti- 
quark pair at infinite separation. Its value is therefore independent of the definition 
of the singlet correlator Gi (r, T) and is related to the renormalized Polyakov loop 



Lren{T) - exp(-^^oo (T)/(2r)) [29]. 

At leading order Fi{r,T) - F^(T) is of Yukawa form (c.f. Eq. ([Ml))- Therefore 
it is useful to define a quantity called the screening function 



S{r,T) = r-{Fi{r,T)-F^{T)). 



(34) 
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Fig. 4. The pre-exponential factor of the color singlet correlators as function of distance r for 
/3 = 2.5 (top) and f3 = 2.7 (bottom). Shown are results for unsmeared spatial links and 10 and 20 
steps of APE smearing. 



This quantity shows exponential decay at distances r > l/T in the deconfined 
phase both in pure gauge theories [18, 32] and full QCD [34]. From its exponential 
decay the Debye screening mass has been determined and turns out to be about 
40% larger than the leading order perturbative value. 

It is interesting to study the screening function using the free energy determined 
from cyclic Wilson loops and make comparison with Coulomb gauge results. This 
analysis has been recently done in Ref. [41] for SU{2) gauge theory. The behavior 
of the screening function at different temperatures is shown in Fig. \5\ At short 
distances (rT < 0.5) the singlet free energy does not depend on the smearing level. 
Furthermore, it is very close to the free energy calculated in Coulomb gauge. At large 
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distances the screening function S{r, T) shows an exponential decay determined 
by a temperature dependent screening mass mi(T), which is equal to the leading 
order Debye mass up to the non-perturbative g"^ corrections: mi — mo + 0{g'^) 
[44, 45]. There is some dependence on the smearing level at larger distances which, 
however, disappears at high temperatures and with increasing the smearing level. 
In particular, for [3 ^ 2.5 it turns out that there is no dependence on smearing level 
for 5 or more smearing steps. For (3 = 2.7 10-20 steps are needed, depending on 
the temperature. Fitting the large distance behavior of the screening function by 
an exponential form exp(— mi(T)r) allows to determine the screening mass mi(T). 
In the inset of Fig. [5] the color singlet screening masses extracted from the fits are 
shown in comparison with the results obtained in Coulomb gauge in Ref. [18]. The 
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solid line is the leading order Debye mass calculated using 2-loop gauge coupling 
= 27rT) in Af5'-scheme. As we see from the figure the screening masses are 
smaller than those calculated in Coulomb gauge and agree well with the leading 
order perturbative prediction. 

4.5. Color adjoint free energy 

The structure of Eqs. (|20p and ([2T|l shows that color adjoint correlator is given 
by a difference of the color averaged and singlet correlators. The color adjoint 
correlator in Coulomb gauge has been studied in pure gauge theory [7, 18], 3-flavor 
QCD [33] and 2-flavor QCD [34]. At low temperatures the color adjoint free energy 
turned out to be significantly smaller than the first hybrid potential contrary to 
the expectations. In fact at sufficiently large distance it was found to be identical 
to the singlet free energy. As has been pointed out in Ref. [19] this is due to the 
non-trivial r-dependence of the overlap factor ci (r) and its deviation from unity. 

We have shown in section [1751 that deviations of the overlap factor ci(r) from 
unity can be greatly reduced when one uses Wilson loops with the smeared spatial 
gauge connection. Therefore it is interesting to see how the adjoint free energy 
behaves in this approach. The numerical analysis has been done in SU{2) gauge 
theory [41] therefore below we will refer to the adjoint free energy as the triplet 
free energy. If we assume that only two states contribute to the Eqs. ((30)) and ((3T|) . 
then from Eq. pT|) it follows that 



with AE{r) — E2{r) — Ei(r). We have seen in section l473l that the temperature 
dependence of the singlet free energy is quite small. In any case it is considerably 
smaller than the temperature dependence of the averaged free energy. Therefore the 
contribution of the excited states to Gi (r, T) is quite small and it is reasonable to 
assume that C2{r) <C 1. From the analysis of the multipole expansion we also expect 
that at small distances, C2(r) ~ (rAgc^))^ [46]. Thus, the temperature dependence 
of F3 (r, T) and its deviation from the hybrid state E2 (r) is due to small deviation 
of Ci(r) from unity. At low temperatures, when AE 3> T these small deviations are 
amplified by the exponential factor. This can be easily verified by subtracting the 
correction Tln(l + i(l — ci)e'^^/"^) from the triplet free energy and assuming that 
El (r) is given by the ground state potential and E2 (f) is given by the first hybrid 
potential as calculated in Ref. [23]. The numerical results are summarized in Fig. 
|6] which shows that after this correction is accounted for in the confined phase the 
triplet free energy at low temperatures agrees reasonably well with the first hybrid 
potential. As temperature increases more excited states contribute. In particular, 
at 0.76Tc the value of the triplet free energy can be accounted for by including the 
next hybrid state [23]. However, at 0.95Tc there are large temperature effects, which 
cannot be explained by including the contribution from only few excited states. 




(35) 
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Fig. 6. The triplet free energy at different temperatures calculated at /3 = 2.5. The filled symbols 
correspond to calculations in Coulomb gauge. Also shown is the first hybrid potential calculated 
in [23]. 



In Fig. [6] we also show the triplet free energy above the deconfinenient tem- 
perature compared to the calculations in Coulomb gauge [18]. It turns out to be 
much smaller than in the confined phase and agrees well with Coulomb gauge re- 
sults. This means that the small deviation of the overlap factor ci (r) from unity 
is unimportant in this case. The triplet free energy monotonically decreases with 
increasing temperature as expected in HTL perturbation theory (c.f. Eq. (|25p ). In 
the limit of high temperatures and short distances, r <C 1/T we have E2{r) = 
as/(4r), AE{r) — as/r,C2{r) ~ Oandci(r) = 1-1-0(0;^). Therefore we can expand 
the logarithm in Eq. (|35|1 to get 

1 fy 

F3{r,T)^+-^ + 0{alT) + 0{a,mD). (36) 
4 r 

Thus the correction due to Ci(r) ^ is much smaller than the expected leading 
order thermal effects in the triplet free energy. 

The strong temperature dependence of the adjoint free energy has been also 
observed in SU{3) gauge theory [7] and in full QCD [33, 34]. Although the ad- 
joint free energy is strongly temperature dependent even at short distances its 
r-dependence is almost the same as r-dependence of the singlet free energy. The 
derivative of the adjoint free energy with respect to r is identical withing errors 
to that of the singlet free energy at short distances up to the Casimir factor, i.e. 
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Table 1. Meson states in different channels 



Fi{r,T)/F[{r,T)^-l/{N^-l) [7]. 

5. Quarkonium spectral functions 

5.1. Meson correlators and spectral functions 

Now we focus the discussion on the relation between the Euclidean meson cor- 
relators and spectral functions at finite temperature. The zero temperature limit 
is straightforward. Most dynamic properties of the finite temperature system are 
incorporated in the spectral function. The spectral function crHiPo-P) for a given 
mesonic channel _ff in a system at temperature T can be defined through the Fourier 
transform of the real time two point functions Z)^ and D'^ or equivalently as the 
imaginary part of the Fourier transformed retarded correlation function [47] , 

Ztt 

= -ImD§{po,p) 

TT 

dT\p,.p) - j ^e^^-D%^<\xo,x) (37) 
D^{xo,x) = (J_ff (xo,x), J//(0,0)) 

D<{xo,x) = {J„{0,0),Jh{xo,x)),xo>0 (38) 

In essence <7h is the Fourier transformation of the thermal average of the commu- 
tator [J{x), J(0)]. 

In the present paper we study local meson operators of the form (c.f. (fT^ with 
U = I) 

Jnit, x) = q{t, x)THqit, x) (39) 
with q a continuous real-time fcrmion position operator and 

= 1 , 75 , 7/. , 757a< 7 7m7i' (40) 

for scalar, pseudo-scalar, vector, axial-vector and tensor channels. The relation of 
these quantum number channels to different meson states is given in Tab. [TJ 
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The correlators D^^'^\xo,x) satisfy the well-known Kubo-Martin-Schwinger 
(KMS) condition [47] 

D>{xo,x)^D<{x„+ilT,x). (41) 
Inserting a complete set of states and using Eq. (|4ip . one gets the expansion 

oh{p,,p) = ¥E™,„(e-^"/^±e-^-"/^) X 

{n\JHmm)?5\p^--kl + k^) (42) 

where Z is the partition function, and fc"^™) refers to the four- moment a of the state 
|n(m)). 

A stable mesonic state contributes a 5 function-like peak to the spectral function: 

^^(Pcp) = mJH\H)\^e{pa)5{p' - ml), (43) 

where niH is the mass of the state and e(po) is the sign function. For a quasi- 
particle in the medium one gets a smeared peak, with the width being the thermal 
width. As one increases the temperature the width increases and at sufficiently 
high temperatures, the contribution from the meson state in the spectral function 
may be sufficiently broad so that it is not very meaningful to speak of it as a well 
defined state any more. The spectral function as defined in Eq. (|42)) can be directly 
accessible by high energy heavy ion experiments. For example, the spectral function 
for the vector current is directly related to the differential thermal cross section for 
the production of dilepton pairs [48] : 

dW 5al^ 1 / ^ /..N 

d^p = 27^pliePofT^if-(P^^P^- (44) 

Then presence or absence of a bound state in the spectral function will manifest 
itself in the peak structure of the differential dilepton rate. 

In finite temperature lattice calculations, one calculates Euclidean time propa- 
gators, usually projected to a given spatial momentum: 



GHir,p)^ J d'xe'P-^TrJH{T,x)JHiO,0)) (45) 

This quantity is an analytical continuation of D^{xo,p) 

Gh{t,p) ^ D>{-iT,p). (46) 

Using this equation and the KMS condition one can easily show that Gh{t,p) is 
related to the spectral function, Eq. (|37p . by an integral equation (see e.g. appendix 
B of Ref. [49]): 



Gh{t,p)= / duja{uj,p)HK{uj,T) 
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This equation is the basic equation for extracting the spectral function from meson 
correlators. Equation (|47p is valid in the continuum. Formally the same spectral 
representation can be written for the Euclidean correlator calculated on the lattice 
G^*(t,p). The corresponding spectral function, however, will be distorted by the 
effect of the finite lattice spacing. These distortions have been calculated in the 
free theory [50, 51]. When discussing the numerical results in following sections the 
subscript H denoting different channels for meson correlators and spectral functions 
will be omitted. 



5.2. Lattice formulations for charmonium physics 

The quarkonium system at zero and finite temperature was studied in [52-55] using 
Wilson-type fermions 



Wilson ^'SW 



tpix) 



where the Dirac operator is defined as 



7-iWilson v7 A 



(48) 



(49) 



with 



Uf,{x)t(j{x + n)- Ul{x - fj.)^p{x - 
Uf,{x)tl){x + ^ji) + Ul{x- - m) - 2V'(a;) 



Furthermore, cr^ ,^ — {7^^,71/} and the field strength tensor is defined as 

4g^,(a:) = U^{x)U,{x + fi)Ul{x + i))Ut{x) + 
U^{x)Ul{x - fl + v)Ul{x - fl)Uf,{x - A) + 
Uj^ix - jji)Ul{x - A - v)U^{x - (1- v)Uu{x - i>) + 
Ul{x - i')U^{x - v)U,.{x + fi- vWlix). 



(50) 



(51) 



(52) 



The last term in the brackets in Eq. ((48|) helps to suppress 0{a) lattice artifacts and 
is called the clover term. Formulations with csw 7^ usually referred to as clover ac- 
tion. The standard Wilson action for fermions corresponds to csw = 0. In the lattice 
literature usually the form with the hopping parameter k — \/{2amQ + 8) is used. 
In this form the fields need to be accordingly normalized ip{x) — > i\)(x)a~'^l'^ j 

In [55] the anisotropic Fermilab formulation is used for the heavy quarks. The 
anisotropy allows to use fine temporal lattice spacings without significant increase 
in the computational cost. The study uses the quenched approximation and the 
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standard Wilson action in the gauge sector for which the relation between the bare 
^0 and the renormalized anisotropy ^ — as/ at is known in a wide range of the gauge 
couphng /3 = 6/5^ [56]. 

The anisotropic clover action [57] is 



Because the lattice spacings in space and time directions are different the spatial 
and temporal Dirac operators as well as the clover terms have different coefficients. 
By tuning the clover coefficients according to the Fermilab prescription [57, 58] it is 
possible to reduce 0{atm()) discretization errors. We will refer to this formulation 
as anisotropic Fermilab formulation. 

In the following we will show results from [55], where anisotropic Fermilab for- 
mulation was used as well as results from the study that used isotropic clover 
action with non-perturbatively determined values of csw [53]. In the study with 
anisotropic Fermilab action two values for the renormalized anisotropy ^ = 2, 4 as 
well as /? = 5.7, 5.9, 6.1, 6.5 were used. These parameters correspond to temporal 
lattice spacings a^^ = 1.905 — 14.12 GeV. To set the scale for the lattice spac- 
ing the traditional phenomenological value rg = 0.5 fm for the Sommer scale [38] 
was used. The Sommer scale rg has also been calculated for anisotropic Wilson 
action for (3 — 5.5 — 6.1 [59]. Alternatively one can estimate the lattice spacing 
from the difference between the mass of ^Pi state and the spin averaged IS" mass: 
AM(^Pi — IS*). To a very good approximation this mass difference is not affected 
by fine and hyperfine splitting and thus is not very sensitive to quenching errors. 
It was found that close to the continuum limit the lattice spacing determined from 
AM{}Pi - IS) is different from that determined from by 10% [57, 60] if the 
phenomenological value r = 0.5 fm is used. Using the value of tq = 0.469(7) de- 
termined in full QCD [61] would give a value for AM(^Pi — IS') splitting which is 
closer to the experimental one, however, the AM(2S'— IS) splitting would be even 
further away from the experimental value [60] . This problem is due to the quenched 
approximation. 

In the studies with isotropic clover action the following values of the lattice 
gauge couphngs have been used (3 = 6.499, 6.640 and 7.192 [53, 62-64]. These 
correspond to lattice spacings a^^ — 4.04, 4.86 and 9.72GeV respectively if the 
value -/a = 420MeV for the string tension is used. 

The continuum meson current in Eq. (|39p . Jh is related to the lattice current 

as 



iWilson 



X ]_ s 




(53) 



Jh Znalil^THtp, 



(54) 



where ip is the lattice quark field in Eq. (|48|) . The renormalization constant Zh 
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can be calculated in perturbation theory or non-perturbatively. For isotropic clover 
action this has been done with both methods (see Ref. [53] and references therein). 

5.3. Bayesian analysis of meson correlators 

The obvious difficulty in the reconstruction of the spectral function from Eq. ([TT)) 
is the fact that the Euclidean correlator is calculated only at 0(10) data points on 
the lattice, while for a reasonable discretization of the integral in Eq. (|47p we need 
C(IOO) degrees of freedom. The problem can be solved using Bayesian analysis of the 
correlator, where one looks for a spectral function which maximizes the conditional 
probability P[a\DH] of having the spectral function a given the data D and some 
prior knowledge H (for reviews see [65, 66]). Different Bayesian methods differ in 
the choice of the prior knowledge. One version of this analysis which is extensively 
used in the literature is the Maximum Entropy Method (MEM) [67, 68] . It has been 
used to study different correlation functions in Quantum Field Theory at zero and 
finite temperature [52-54, 62, 65, 68-79]. In this method the basic prior knowledge 
is the positivity of the spectral function and the prior knowledge is given by the 
Shannon - Janes entropy 



The real function m(uj) is called the default model and parametrizes all additional 
prior knowledge about the spectral functions, e.g. such as the asymptotic behavior 
at high energy [65, 68]. For this case the conditional probability becomes 



with being the standard likelihood function and a a real parameter. Previously 
in the MEM analysis of the meson spectral functions the Bryan's algorithm was 
used [67] . A new algorithm was introduced in [55] . It is worth to make connection 
between this method and the Bryan algorithm. In both cases the true problem is 
number-of-data dimensional - in more dimensions the problem would be under- 
determined. To find the relevant subspace, the Bryan algorithm uses singular value 
decomposition, while the new algorithm finds the same relevant subspace by exact 
mathematical transformations. Although the method of identifying the subspace is 
different, the result is the same, and in both cases one proceeds with solving the 
original problem in this restricted subspace. The advantage of the new algorithm 
is that it is more stable numerically when one reconstructs quarkonium spectral 
functions at zero temperature. 

The comparison of the two algorithms was done for the pseudo-scalar spectral 
function for (3 = 6.1, ^ = 4 [55]. The problem with the Bryan algorithm is that 
it does not work well for charmonium correlators if the time extent is sufficiently 
large, which is the case at low temperatures; the iterative procedure does not always 
converge. For instance at (3 — 6.1, ^ = 4 and 16'^ x 96 lattice it was possible to get 




P[a\DH]=exp{--x^ + aS), 



(55) 
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Fig. 7. The pseudo-scalar spectral function at zero temperature for three finest lattice spacings. 
The horizontal line corresponds to the spectral function in the free massless limit at zero lattice 
spacing. 



the spectral functions using the Bryan algorithm only when using Tmax = 24 data 
points in the time direction. With the new algorithm there is no restriction on Tmax 
which can be as large as Nt/2. 



5.4. Charmonium spectral functions at zero temperature 

5.4.1. Pseudo-scalar and vector spectral functions at zero temperature 

In this subsection we discuss the results of Ref. [55] on charmonium spectral func- 
tions obtained using MEM. The zero temperature spectral functions for three dif- 
ferent lattice spacings obtained by the new method from Ref. [55] are summarized 
in Fig. [T] Here the simple default model m{uj) — 1 (in units of the spatial lattice 
spacing) is used. To get a feeling for the statistical errors in the spectral functions 
its mean value in some interval / is calculated: 

- Ii dujgjuj) 

Then the error on a is calculated using standard jackknife method. These errors are 
shown in Fig. [71 where the length of the intervals are shown as horizontal error bars. 
As one can see from the figure, the ric{lS) can be identified very well. The second 
peak is likely to correspond to excited states. Because of the heavy quark mass 
the splitting between different radial excitations is small and MEM cannot resolve 
different excitations individually but rather produces a second broad peak to which 
all radial excitation contribute. This can be seen from the fact that the amplitude 
of the second peak (i.e. the area under the peak) is more than two times larger than 
the first one. Physical considerations tell us that it should be smaller than the first 
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Fig. 8. The default model dependence of the pseudo-scalar spectral function at the finest lattice 
spacing (/3 = 6.5). In the inset the default models corresponding to different spectral functions 
are shown. 



amplitude if it was a 2S state. When comparing amplitudes and peak positions from 
MEM analysis and from double exponential fits a very good agreement for the first 
peak and a fair agreement for the second peak are found. This gives confidence that 
at zero temperature charmonium properties can be reproduced well with MEM. 

For energies larger than 5GeV one probably sees a continuum in the spectral 
functions which is distorted by finite lattice spacing. In particular the spectral 
function is zero above some energy which scales roughly as aj^. Note that for 
uj < 5GeV the spectral function does not depend on the lattice spacing. 

One should control how the result depends on the default model. In Fig. [5] we 
show the spectral function for three different default models. One can see that the 
default model dependence is significant only for w > 5 GeV. This is not surprising 
as there are very few time slices which are sensitive to the spectral functions at 
a; > 5 GeV, while the first peak is well determined by the large distance behavior 
of the correlator. 

The spectral function in the vector channel defined as 

crv(w) = ^^crjj(w), (57) 

i 

was also calculated and is shown in Fig. [9] for the three finest lattice spacings. The 
conclusions which can be derived from this figure are similar to the ones discussed 
above for the pseudo-scalar channel. The first peak corresponds to the J/il){lS) 
state, the second peak most likely is a combination of 2S and higher excited states, 
finally there is a continuum above 5 GeV which is, however, distorted by lattice 
artifacts. The similarity between the pseudo-scalar and vector channel is, of course, 
expected. The lower lying states in theses channels differ only by small hyperfine 
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Fig. 9. The vector spectral function at zero temperature for three finest lattice spacings. The 
horizontal line is the spectral function in the free massless limit. 



splitting. 



5.4.2. Spectral functions for P states 

The spectral functions in the scalar, axial-vector and tensor channels which have 
the IP charmonia as the ground state were also calculated in Ref. [55]. The scalar 
spectral functions reconstructed using MEM are shown in Fig. [TOl The first peak 
corresponds to XcO state, but it is not resolved as well as the ground state in the 
pseudo-scalar channel. This is due to the fact that the scalar correlator is consider- 
ably more noisy than the pseudo-scalar or vector correlator. This can be understood 
as follows. For the heavy quark mass the contribution of the ground state in the 
scalar channel is suppressed as 1/m^ relative to the ground state contribution in 
the pseudo-scalar and vector channels, and therefore it is considerably smaller than 
the continuum contribution to the scalar correlator. For the two finest lattice spac- 
ings there is a second peak which may correspond to a combination of excited P 
states. Above w > 5 GeV we see a continuum which is strongly distorted by lattice 
artifacts and probably also by MEM. 

The spectral functions in the axial- vector and tensor channels are shown in Fig. 
im They look similar to the scalar spectral functions. As in the scalar channel 
the first peak is less pronounced than in the case of S-wa.ve charmonium spectral 
functions, and it corresponds to Xci and he state, respectively. The continuum part 
of the spectral function is again strongly distorted. 
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Fig. 10. The scalar spectral function at zero temperature for three finest lattice spacings. 
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Fig. 11. The axial-vector (left) and tensor (right) spectral functions at zero temperature for three 
lattice spacings. 



5.5. Charmonium correlators at finite temperature 

The reconstruction of the spectral functions from lattice correlators is difficult al- 
ready at zero temperatures. At finite temperature it is even more difficult to control 
the systematic errors in the spectral functions reconstructed from MEM. This is 
because with increasing temperature the maximal time extent Tmax is decreasing 
as 1/T. Also the number of data points available for the analysis becomes smaller. 
Therefore other methods which can give some information about the change of the 
spectral functions as the temperature is increasing are often used. The temperature 
dependence of the spectral function will manifest itself in the temperature depen- 
dence of the lattice correlator G(r, T). Looking at Eq. ([47|) it is easy to see that the 
temperature dependence of G(r, T) comes from the temperature dependence of the 
spectral function (t(w, T) and the temperature dependence of the kernel K(t, lu, T). 
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Fig. 12. The ratio GjGr 



for tile pseudo-scalar channei for tiie two finer ^ = 4 lattices. 



To separate out the trivial temperature dependence due to K{t,u}, T) one calculates 
the reconstructed correlator [53] 

/>OC 

G™(r,T)= / dLua{u,T ^Q)K{t,u,T). (58) 
Jo 

If the spectral function does not change with increasing temperature we expect 
G{T,T)/Grecon{T,T) = 1. In Rcf. [55] an extensive study of the temperature de- 
pendence of this ratio for different channels at different lattice spacings was carried 
out. 



5.5.1. The pseudo-scalar correlators 

First we present results from Ref. [55] for the temperature dependence of the pseudo- 
scalar correlators. In Fig.[l2]we show numerical results for G/Grecon on lattices with 
^ = 4. We see almost no change in the pseudo-scalar correlator till temperatures as 
high as 1.2Tt,. The temperature dependence of the pseudo-scalar correlator remains 
small for temperatures below 1.5Tc. Medium modifications of the correlator slowly 
turn on as we increase the temperature above this value. From the figures it is clear 
that the temperature dependence of the correlators is not affected significantly by 
the finite lattice spacing. The very small temperature dependence of the pseudo- 
scalar correlator suggests that the corresponding ground state ric{lS) may survive 
till temperatures as high a 1.5Tc. The temperature dependence of the correlator 
found in this study is similar to that of Ref. [53], where isotropic lattices with 
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very small lattice spacings, a^^ — 4.86,9.72 GeV have been used. However, at 
temperatures higher than 1.5Tc the deviations of G/Grecon from unity become 
slightly larger than those found in Ref. [53]. This is possibly due to the fact that 
cutoff effects are more important at higher temperatures. Thus despite similarities 
of the temperature dependence of the pseudo-scalar correlator to findings of Rcf. [53] 
there are quantitative differences. One should note, however, statistical errors and 
systematic uncertainties are larger in the analysis presented in Ref. [53] than in 
Ref. [55]. In the study on isotropic lattices the ratio G/Grecon starts to depend 
more strongly on the temperature only around 3Tc [53] . 

5.5.2. The P-wave correlators 

Next we present the temperature dependence of the scalar, axial-vector and tensor 
correlators corresponding to P-states. The r dependence of the scalar correlator 
was studied on ^ = 2,4 lattices [55]. The numerical results on fine lattices with 
^ = 4 are shown in Fig. [121 We see some differences in G/Grecon calculated at 
P = 6.1 and f3 — 6.5. Thus the cutoff dependence of G/Grecon is larger in the scalar 
channel than in the pseudo-scalar one. For /3 = 6.1 and ^ = 4 the calculations were 
done on 24^^ x 24 lattice to check finite volume effects. The corresponding results 
are shown in Fig. 1131 indicating that the finite volume effects are small. On the 
finest lattice the enhancement of the scalar correlator is very similar to that found 
in calculations done on isotropic lattices [53] , but small quantitative differences can 
be identified. 

In Figs . [Ml and [TSl we show the temperature dependence of the axial- vector and 
tensor correlators respectively for ^ = 4. Qualitatively their behavior is very similar 
to the scalar correlator but the enhancement over the zero temperature result is 
larger. The results for the axial-vector correlators again are very similar to those 
published in Ref. [53]. The difference in G/Grecon calculated at /3 = 6.1 and f3 = 6.5 
are smaller than in the scalar channel. 

The large increase in the scalar, axial-vector and tensor correlators may be in- 
terpreted as indicator of strong modification of the corresponding spectral function 
and, possibly the dissolution of IP charmonia states. However, as we will see in the 
next sections the situation is more complicated. It has been noticed in [80] that the 
increase in G/Grecon may be due to the zero mode contribution. 

5.5.3. The vector correlator 

The numerical results for the vector correlator of Ref. [55] are shown in Fig. [16] 
for ^ = 4. As one can see from the figures the temperature dependence of G/Grec 
is different from the pseudo-scalar case and this ratio is larger than unity for all 
lattice spacings. Similar results have been obtained on isotropic lattices [63]. The 
enhancement of the vector correlator is due to the presence of the transport contri- 
bution in the spectral function [49, 81] and will be discussed in section (5. 101 Since 
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Fig. 14. The ratio G/Gr econ for the axial-vector channel for ^ — 4 lattices. 
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Fig. 17. The dependence of the reconstructed pseudo-scalar spectral function on the maximal 
temporal extent for f) = 6.5. In the analysis the default model m(a;) = 1 has been used. 

the vector current is conserved the vector correlator also caries information about 
the transport of heavy quarks in the plasma. This shows up as a peek in the spectral 
function at zero energy, leading to the observed enhancement in G / Grecon- 

5.6. Charmonium spectral functions at finite temperature 

In section 15.41 we have seen that using MEM one can reconstruct well the main 
features of the spectral function, in particular the ground state properties. At finite 
temperature the situation becomes worse because the temporal extent is decreasing. 
The maximal time separation is Tmax = 1/(27"). As a consequence it is no longer 
possible to isolate the ground state well. Also the number of available data points 
becomes smaller. While the later limitation can be overcome by using smaller and 
smaller lattice spacings in time direction the former limitation is always present. 
Therefore we should investigate systematic effects due to limited extent of the 
temporal direction. It appears that the pseudo-scalar channel is the most suitable 
case for this investigation as at zero temperature it is well under control and there 
is no contribution from heavy quark transport. To estimate the effect of limited 
temporal extent in Ref. [55] the spectral functions at zero temperature is calculated 
considering only the first Ndata time-slices in the analysis for j3 = 6.5, ^ = 4. The 
result of this calculation is shown in Fig. [17] where Ndata = 80, 40, 20 and 16. 
The last two values correspond to the finite temperature lattices in the deconfined 
phase. In this case we see the first peak quite clearly. As one can see from the figure 
already for Ndata — 40 and T,nax = 0.56 fm the second peak corresponding to radial 
excitation is no longer visible and the first peak becomes significantly broader. The 
position of the first peak, however, is unchanged. As the number of data points is 
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further decreased {Ndata — 20, 16) we see further broadening of the first peak and a 
smah shift of the peak position to higher energies. These systematic effects should 
be taken into account when analyzing the spectral functions at finite temperature. 
Therefore when studying the spectral functions at finite temperature we always 
compare with the zero temperature spectral functions reconstructed with the same 
number of data points and Tmax as available at that temperature. 

In Fig. [18] the spectral functions at different temperatures are shown together 
with the zero temperature spectral functions [55]. As a default model for T = 1.2rc 
and T — 1.5Tc m{uj) — 0.01 is used. For T — 2.0Tc m{uj) = 1 is used since the use of 
m{uj) — 0.01 resulted in numerical problems in the MEM analysis. The figure shows 
that the pseudo-scalar spectral function is not modified till 1.5Tc within the errors 
of the calculations. This is consistent with the conclusions of Ref. [53, 54]. One 
should note, however, that it is difficult to make any conclusive statement based on 
the shape of the spectral functions as this was done in the above mentioned works. 
The dependence of the reconstructed spectral functions on the default model 'm{uj) 
is much stronger at finite temperature. The spectral functions were reconstructed 
using different types of default models. For all temperatures T < 1.5Tc the difference 
between the finite temperature spectral function and the zero temperature one is 
very small compared to the statistical errors for all default model considered here. 
In particular a use is made of the default models constructed from the high energy 
part of the lattice spectral functions calculated at zero temperature as this was 
done in Ref. [53] . The idea is that at sufficiently high energy the spectral function 
is dominated by the continuum and is temperature independent. Therefore it is 
suitable to provide the prior knowledge, i.e. the default model. With this default 
model the spectral functions were calculated at T = (1.07 — 1.5)rc [55]. Very little 
temperature dependence of the spectral functions was found, see Fig. [111 Note, 
however, that for this choice of the default model no clear peak can be identified in 
the spectral functions if T > 1.2rc- 

The spectral function in the vector channel is also calculated [55] . The results are 
shown in Fig. [20] for the default model m{uj) — 0.01 . As this was already discussed 
in the previous section the basic difference between the pseudo-scalar and vector 
spectral functions at finite temperature is the presence of the transport peak at 
w ~ 0. The difference of the temperature dependence of the vector and pseudo- 
scalar correlators is consistent with this assumption. The vector spectral function 
reconstructed with MEM shows no evidence of the transport peak at a; ~ 0. On the 
other hand the spectral function at 1.2rc differs from the zero temperature spectral 
function, in particular the first peak is shifted to smaller co values. We believe that 
this is a problem of the MEM analysis which cannot resolve the peak at a; ~ but 
instead mimics its effect by shifting the J/ip peak to smaller uj. Also at 2.4Tc the 
spectral function extends to smaller u values than in the pseudo-scalar correlator 
which again indicates some structure at a; ~ 0. The analysis of the vector spectral 
functions using other choices for the default model always indicates that the spectral 
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Fig. 18. The pseudo-scalar spectral function for f} = 6.5 and Nt = 40, 32, 20 corresponding to 
temperatures 1.2Tc, 1.5Tc and 2.0rc. In the analysis the default model m(a;) = 0.01 has been 
used for 1.2Tc and 1.5Tc, while for 2.0rc we used m(u)) = 1. 
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Fig. 19. The pseudo-scalar spectral function at different temperatures together with the zero 
temperature spectral functions reconstructed using default model coming from the high energy 
part of the zero temperature spectral function. 



functions at finite temperature differs from the zero temperature spectral functions 
and extend to significantly smaller lo values. 



5.7. Charmonium correlators and spectral functions at finite 
momenta 

So far we reviewed charmonia at zero spatial momentum, i.e. charmonia at rest 
in the heatbath's rest frame. It is certainly of interest to study the temperature 
dependence of correlators and spectral functions at non-zero spatial momentum. 
Such a study has been done using isotropic lattices with lattice spacing = 
4.86GeV and 9.72 GeV [63]. It has been found that the pseudo-scalar correlators are 
enhanced compared to the zero temperature correlators for non-vanishing spatial 
momenta, see Fig. 1211 Furthermore, the enhancement of vector correlator at finite 
spatial momentum is larger than at zero spatial momentum. 

In Ref. [55] the finite momentum pseudo-scalar correlators are calculated on 
anisotropic lattices at /3 = 6.1, ^ = 4 for different temperatures. The differences 
in G/Grecon Calculated in this work and in Refs. [53, 63] are present already at 
zero momentum and are presumably due to finite lattice spacing errors. Apart 
from this the momentum dependence of the pseudo-scalar correlators is similar to 
the findings of Refs. [53, 63]. It would be interesting to see if the difference in 
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Fig. 20. The vector spectral function for /3 = 6.5 and Nt = 40, 20 corresponding to temperatures 
1.2rc and 2.4rc. In the analysis the default model m{u}) = 0.01 has been used. 



the temperature dependence of the correlators at zero and finite spatial momenta 
is due to a contribution to the spectral functions below the light cone at finite 
temperature [50, 51]. 

5.8. Bottomonium spectral functions at zero temperature 

The use of Fermilab formulation described in the previous sections allows for a 
study of bottomonium for the same range of lattice spacings. Usually bottomonium 
is studied using lattice NRQCD (see e.g Ref. [82]). First study of bottomonium 
within the relativistic framework was presented in Ref. [83]. More recently it was 
studied in [55], where as before the lattice spacing is fixed by the Sommer scale ro, 
assuming tq — 0.5fm. Note that the lattice spacings determined from tq are about 
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Fig. 21. Comparison of correlators at 1.1 Tc and 1.5 Tc for pseudoscalar and (transverse) vector 
charmonia with those reconstructed from spectral function at 0.75 Tc. 



20% smaller than in Ref. [83] where it was determined from bottomonium ^ Pi — 15 
mass splitting. As the result the new estimates of the T mass are smaller than those 
in Ref. [83]. However, one finds good agreement if the values of the lattice spacing 
quoted in Ref. [83] are used to calculate the physical masses. 

Using MEM the spectral functions in different channels for three lattice spac- 
ings were analyzed. In Fig. [22] we show the spectral functions in the pseudo-scalar 
channel. Since the physical quark mass is different at different lattice spacings the 
horizontal scale was shifted by the difference of the calculated T-mass and the 
corresponding experimental value. We can see that the first peak in the spectral 
function corresponds to the ?7f,(lS') state and its position is independent of the lat- 
tice spacing. The remaining details of the spectral functions are cut-off dependent 
and we cannot distinguish the excited states from the continuum. The position and 
the amplitude of the first peak in the spectral functions is in good agreement with 
the results of simple exponential fit. As in the charmonium case the maximal en- 
ergy Umax for which the spectral function is non-zero scales approximately as aj^. 
Similar results have been obtained in the vector channel. 

The spectral function in the scalar channel is shown in Fig. [23l As it was the 
case for charmonium the correlators in this channel are more noisy than in the 
pseudo-scalar channel and as a result it is more difficult to reconstruct the spectral 
function. Nevertheless we are able to reconstruct the Xfco state which is the first 
peak in the spectral function. The peak position and the amplitude are in reasonable 
agreement with the result of the simple exponential fit. 

5.9. Bottomonium at finite temperature 

Besides calculating the bottomonium spectral function at zero temperature Ref. [55] 
presents a study of the temperature dependence of bottomonium correlators to see 
medium modification of bottomonia properties. In Fig. [24] we show G/Grecon for 
vector and pseudo-scalar channel at different lattice spacings. This ratio appears 
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Fig. 22. The pseudo-scalar bottomonium spectral function at zero temperature for different lattice 
spacings. 




to be temperature independent and very close to unity up to quite high tempera- 
tures. This is consistent with the expectation that IS bottomonia are smaller than 
IS charmonia and thus are less effected by the medium. They could survive till 
significantly higher temperatures. Compared to charmonium case the difference be- 
tween the pseudo-scalar and vector channels is smaller. This is also expected as 
the transport contribution which is responsible for this difference is proportional to 
~ exp(— mc,{,/r), and thus is much smaller for bottom quarks (see the discussion in 
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in the pseudo-scalar (top) and vector channels at different lattice 



the next section). Similar temperature dependence of the pseudo-scalar bottomo- 
nium correlator has been found in calculations with isotropic clover action [84]. 

The temperature dependence of the scalar correlator is shown in Fig. [^Sl Con- 
trary to the pseudo-scalar and vector correlators it shows strong temperature de- 
pendence and G/Grecon is significantly larger than unity already at l.lTc. Again, 
similar enhancement in G/Grecon has been observed in isotropic lattice calcula- 
tions [84]. 

The spectral functions were reconstructed at finite temperature for /3 = 6.3. For 
the pseudo-scalar channel the results are shown in Fig. [521 As in the charmonium 
case we compare the finite temperature spectral function with the zero temperature 
spectral function obtained with the same number of data points and time interval. 
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Fig. 26. The pseudo-scalar bottomonium spectral function at finite temperature. 

As expected the spectral function shows no temperature dependence within errors. 
On the other hand it was not possible to reliably reconstruct the scalar spectral 
function at finite temperature due to numerical problems. Presumably much more 
statistics is needed to get some information about the scalar spectral function. 

5.10. Zero modes contribution 

At finite temperature quarkonium spectral functions contain information about 
states containing a quark anti-quark pair as well as scattering of the external probe 
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Fig. 27. The ratio of the derivatives G' {t,T) /G'^^^^„{t,T) in the scalar channel (left) and axial- 
vector channel (right) calculated at /3 = 7.192. The results from anisotropic lattice calculations at 
/3 = 6.5 [55] are also shown (open symbols). 



off a heavy quark from the medium. The later gives a contribution to the spectral 
function below the light cone (w < k). In the limit of zero momentum it becomes 
X^{T)uj6{uj) in the free theory. The generalized susceptibilities X^{T) were calculated 
in Ref. [51] in the free theory. Interaction with the medium leads to the broadening 
of the delta function, which becomes a Lorentzian with a small width [81]. Because 
the quark anti-quark pair contributes to the spectral function at energies co > 2m 
it is reasonable to separate quarkonium spectral function into two terms [64] 

a\u;,T) = al,^^{u;,T) + al^{cu,T). (59) 

Here CThighC'^j 2^) is the high energy part of the spectral functions which is non- 
zero only for uj > 2m and describes the propagation of bound or unbound quark 
anti-quark pairs. On the other hand al^^{uj,T) receives the dominant contribution 
at a; ~ 0. Because the width of the peak at cj ~ is small the later gives an 
almost constant contribution to the Euclidean correlator, which is called the zero 
mode contribution. We can write an analogous decomposition for the Euclidean 
correlator 

G^T,T) = GLigh(r,r) + GU(r,r). (60) 

To a very good approximation GJq^(t, T) — x^{T)T, i.e. constant. 

In the previous sections we used the ratio G{t,T) / Grecon{T,T) to study the 
temperature dependence of the correlators. This dependence comes separately from 
the high energy part and low energy part, which gives the zero mode contribution. 
The zero mode contribution is absent in the derivative of the correlator with respect 
to T. Therefore one can study the temperature dependence of the correlators induced 
by change of bound state properties and/or its dissolution by considering the ratio 
of the derivatives of the correlators G'(r, T)/GJ,gj,o„(T, T). 

The temperature dependence of scalar (fT3|) and axial-vector correlators (|14|) 
has been presented in previous sections in terms of G/Grecon- It is temperature 
independent in the confined phase and shows large enhancement in the deconfined 
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Fig. 28. The thermal mass correction as function of the temperature. The dashed Hne and the 
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Fig. 29. The thermal velocity of heavy quarks (left) and the zero mode contribution to the axial- 
vector correlator Gf^^/T^ = 'x°'^ jT"^ (right). The lines show the prediction of the quasi-particle 
model with mgg(T). The open symbols show the thermal velocity squared estimated on anisotropic 
lattice s [55]. 



phase. This large enhancement is present both in charmonium and bottomonium 
correlators. To eliminate the zero mode contribution the derivative of the correla- 
tors and the corresponding ratio G"(r, T)/G'J,gco„(T, T) have been calculated in [64] 
using isotropic lattices. The numerical results for this ratio at /3 = 7.192 are shown 
in Fig. [27l The results from anisotropic lattices [55] are also shown. There is good 
agreement between the results obtained from isotropic and anisotropic lattices. As 
one can see from Fig. [27] G"(T,T)/G'^gco„(r,T) shows very little temperature de- 
pendence and is close to unity. This means that almost the entire temperature 
dependence of the scalar and axial-vector correlators is due to zero mode contri- 
bution and G'high(''"j 2^) is temperature independent. The temperature dependence 
of the S-states seen in previous sections is also greatly reduced and the agreement 
between isotropic and anisotropic calculations is better for the ratio of derivatives. 
For the pseudo-scalar channel this behavior can be explained by the presence of a 
small negative zero mode [64]. 
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Since the high energy part of the quarkonium correlators turns out to be temper- 
ature independent to very good approximation one can assume that GJ^jgi^ (r, T) ~ 
G5,ggo„(T, T). Then the low energy part of the correlators, i.e. the zero mode con- 
tribution can be evaluated as Gf„^(r,T) = G*(t,T) - GI.^^^^{t,T) [64]. Current 
lattice data are not precise enough to see a clear r dependence in GJ^^ (r, T) and 
it is compatible with being a constant. Therefore, one could take the value of low 
energy part of the correlator at the midpoint as an estimate for the zero mode 
contribution, GIJt = 1/{2T),T) = Tx^T). 

The zero mode contribution is related to the propagation of single (unbound) 
heavy quark in the medium. Therefore, it is natural to describe its temperature de- 
pendence in terms of a quasi-particle model with effective temperature dependent 
heavy quark masses. Since the temporal component of the vector correlator has no 
high energy part, i.e. G'«=°(r,T) = -Tx"''"(T) = -Tx{T) it is most suitable for 
fixing the effective quark mass. Matching the lattice data to the free theory expres- 
sion for x(T) one can determine the effective quasiparticle masses rricff [64]. The 
results of this analysis are shown in Fig. [28] for different values of the constituent 
heavy quark mass, including the bottom quarks, in terms of thermal mass correction 
5mcs{T) = mcs{T) — m. As one can see from the figure the thermal mass correction 
decreases monotonically with increasing the constituent quark mass and increas- 
ing temperature for bottom quarks it is not incompatible with the leading-order 
perturbative prediction 

dTOcff = -o— ; "i-D, (61) 

o 47r 

with rriD ~ g{T)T being the perturbative Debye mass {Nf = because we work in 
the quenched approximation). 

Having determined the effective heavy quark mass we can study the zero mode 
contribution in other channels. If the quasi-particle model is correct the zero mode 
contribution should be a function of m^ff/T only. Therefore in Fig. [^Hl the tempera- 
ture dependence of zero mode contribution for the vector and axial- vector channel 
is shown as function of T/rrieg. Indeed, all the lattice data seem to fall on one 
curve within errors, which agrees with the quasi-particle model prediction shown 
as the black lines. For the vector channel we show the data in terms of the ratio 
GJ'q^(T)/G"'^°, where GJ'q^(T) is the sum over all spatial components. This is be- 
cause this quantity does not depend on the renormalization [64] and has simple 
physical interpretation in terms of averaged thermal velocity squared u^jj. The later 
follows from the fact that due to the large quark mass the Boltzmann approximation 
can be used and we have 



The zero mode contribution in the vector channel has also been studied on 
anisotropic lattices [55] and in Fig. [29] we also show the corresponding results for 
thermal velocity squared. 
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6. Potential models at finite temperature 

Quarkonium properties at finite temperature have been studied in potential models 
since the famous paper of Matsui and Satz [1] (for recent review see [2]). The basic 
idea behind this is that color screening will modify the potential, which becomes 
short range and cannot support bound states of heavy quarks at sufficiently high 
temperatures. As discussed in section [2] this approach can be justified if there is a 
separation of the scales related to binding energy and other scales in the problem, 
like the inverse size of the bound states and the temperature. Close to the QCD 
transition this separation of scales is not obvious, however, the lattice calculation of 
static quark correlators, discussed in section |4] show that screening effects are very 
strong already in the transition region. In fact, correlation functions of static quark 
anti-quark pair show significant temperature modifications already at distances 
similar to quarkonium size. Therefore, we may expect the most of quarkonium 
bound states dissolve in the deconfined phase at temperatures close to the transition 
temperature. This seemingly contradicts to the small temperature dependence of 
quarkonium correlators and spectral functions discussed in the previous section. 

Due to the heavy quark mass quarkonium spectral functions can be calculated 
in the potential approach by relating the spectral functions in the threshold region 
to the non-relativistic Green function [14, 85, 86] 

a{Lo) = if-ImG"'-(r, P, E)\,_p_^ , (63) 

<j{uj) = i^-ijImV • V'G"'-(f, r', E%^p^, , (64) 

for S'-wave, and P-wave quarkonia, respectively. Here E = u — 2m. The pre-factor 
K accounts for relativistic and radiative corrections. The non-relativistic Green 
function is calculated from the Schrodinger equation with a delta-function on the 
right hand side. Away from the threshold the spectral function is matched to the 
perturbative result. It has been shown that this approach can provide a fair de- 
scription of quarkonium correlators at zero temperature [87]. Other approaches 
to calculate quarkonium spectral functions in potential models were proposed in 
Refs. [49, 88-91]. Close to the transition temperature the perturbative calculations 
of the potential are not applicable. Therefore in Ref. [85] the singlet free energy 
calculated in quenched lattice QCD has been used to construct the potential. The 
drawback of this approach is that there is no one to one correspondence between 
the singlet free energy and the potential in the non-perturbative domain. Therefore 
in Ref. [85] several possible forms of the potential compatible with the lattice data 
have been considered. Furthermore, lattice calculations do not have much informa- 
tion about the imaginary part of the potential. Therefore the imaginary part of 
the potential has been approximated by a small constant term. The results of the 
calculations of the charmonium and bottomonium spectral functions for S'-wave are 
shown in Figl30l In the case of charmonium all bound states are melted at tem- 
peratures higher than l.2Tc. We see, however a significant threshold enhancement. 
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Fig. 30. The charmonium (top) and bottomonium (bottom) spectral functions at different tem- 
peratures. For charmonium we also show the spectral functions from lattice QCD obtained from 
the MEM at l.bTc- The error-bars on the lattice spectral function correspond to the statisti- 
cal error of the spectral function integrated in the oj-interval corresponding to the horizontal 
error-bars. The insets show the corresponding ratio G/Grecon together with the results from 
anisotropic lattice calculations [55]. For charmonium, lattice calculations of G/Grecon are shown 
for T = 1.2Tc (squares), 1.5rc (circles), and 2.0Tc (triangles). For bottomonium lattice data are 
shown for T = 1.5Tc (circles) and l.STc (triangles). 



i.e. near the threshold the spectral function is much larger than in the free case. In 
the bottomonium case only the ground state survives in the deconfined phase. At 
temperatures above 2Tj, we see the melting of the ground state as well. 

Quarkonium spectral functions calculated using the perturbative potential, 
where the imaginary part is fully taken into account, show similar qualitative fea- 
tures [14, 86]. Here too a significant threshold enhancement is seen. The imaginary 
part of the potential plays an important role in weakening the bound state peak 
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or transforming it to mere threshold enhancement. Let us note that the next-to- 
leading order perturbative correction to the quarkonium spectral functions also give 
rise to significant enhancement in the threshold region [92]. 

From the spectral functions we can calculate the quarkonium correlation func- 
tions in Euclidean time G(r, T) and compare them to the available lattice data. This 
comparison is shown in Fig. [30] for the ratio G(r, T)/Grecon{T, T). As one can see 
from the figure the melting of bound states does not lead to large change in the Eu- 
clidean correlation functions. The ratio G{t,T) /GreconiT,T) calculated in the po- 
tential model is flat and temperature independent in agreement with lattice calcula- 
tions. This means that threshold enhancement can compensate for melting of bound 
states in terms of the Euclidean correlators. Spectral functions of P-wave quarko- 
nium have been also calculated [14, 85] and show significant threshold enhancement 
as well. As the consequence the ratio of the derivatives G'(r, T)/ G'^.^^^^{t, T) is tem- 
perature independent and close to unity [85] in agreement with lattice calculations 
shown in the previous section. 

The analysis discussed above has been done in quenched QCD. This is because 
only in quenched QCD we have sufficiently precise lattice calculations of quarko- 
nium correlators. Potential model calculations of the spectral functions have been 
extended to 2-1-1 flavor QCD using the lattice data discussed in section 21 These 
calculations show that all quarkonium states except the ground state bottomo- 
nium dissolve in the quark gluon plasma [93] . The upper limits on the dissociation 
temperatures for different quarkonium states obtained in this analysis are given in 
Tabled 



state 


Xc 


^' 


J/i, 




Xb 


T 










1.2T, 


1.3T, 





Table 2. Upper bound on the dissociation temperatures for different quarkonium states in 2+1 
flavor QCD [93]. 



7. Conclusion 

In this paper we discussed quarkonium properties in quark gluon plasma. We dis- 
cussed how color screening can be studied non-perturbatively on the lattice using 
spatial correlation functions of static quark and anti-quark. We discussed the im- 
portance of singlet and adjoint (triplet for SU{2) or octet for SU{3)) degrees of 
freedom for understanding the temperature dependence of static meson (quark- 
antiquark) correlators. It has been shown how color singlet and adjoint degrees of 
freedom can be defined in the effective field theory framework, the so-called thermal 
pNRQCD. 

We have seen that the singlet correlators of static quark and anti-quark show 
strong screening effects at distances comparable to quarkonium size. Thus it is nat- 
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ural to expect that most quarkonium states melt in quark gluon plasma. Quite 
surprisingly lattice calculations of quarkonium correlators in Euclidean time show 
very little temperature dependence. Wc reviewed the current status of these calcu- 
lations in quenched QCD. It turns out that the spatial lattice spacing should be 
smaller than (4GeV)~^ to have full control over the discretization errors. Therefore 
the extension of these calculations to full QCD pioneered in Ref. [94] will be quite 
difficult. 

The only source of significant temperature dependence of quarkonium corre- 
lators is the zero mode contribution which is not related to bound states but to 
transport properties of heavy quarks in the medium [81]. We have shown that 
this contribution can be well described by a quasi-particle model with temperature 
dependent heavy quark mass. We have presented the calculations of quarkonium 
spectral functions in potential model and have shown that all quarkonium states, 
except the ground state bottomonium melt in the dcconfined phase. We have also 
shown how the seemingly existing contradiction of strong color screening leading 
to quarkonium melting and very weak temperature dependence of quarkonium cor- 
relators can be resolved within potential models. It turns out that strong thresh- 
old enhancement of quarkonium spectral functions can compensate the absence 
of bound state and result in Euclidean correlarors, which are almost temperature 
independent. 

The fact that no charmonium bound states can exist in quark gluon plasma 
has important consequences for describing charmonium production in heavy ion 
collisions. In absence of bound states heavy quarks in the plasma can be treated 
quasi-classically using Langevin dynamics. In this scenario the residual correlation 
between the heavy quark and anti-quark, visible in the threshold enhancement and 
the finite life time of the plasma play an important role [95] . It turns out that using 
these ideas and the interactions of heavy quark determined in lattice QCD it is 
possible to uderstand the J/V' suppression pattern at RHIC [95] . 
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